Controlling integrability in a quasi- ID atom-dimer mixture 
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We analytically study the atom-dimer scattering problem in the near-integrable limit when the 
oscillator length Iq of the transverse confinement is smaller than the dimer size, ~ ^o/loili where a < 
is the interatomic scattering length. The leading contributions to the atom-diatom reflection and 
break-up probabilities are proportional to a in the bosonic case and to a for the t-ti scattering in 
a two-component fermionic mixture. We show that by tuning a and Iq one can control the "degree of 
integrability" in a quasi-lD atom-dimer mixture in an extremely wide range leaving thermodynamic 
quantities unchanged. We find that the relaxation to deeply bound states in the fermionic (bosonic) 
case is slower (faster) than transitions between different Bethe ansatz states. We propose a realistic 
experiment for detailed studies of the crossover from integrable to nonintegrable dynamics. 

PACS numbers: 34.50.-s, 67.85. Pq 



I. INTRODUCTION 

Ultracold gases allow unprecedented experimental con- 
trol over key parameters of a many-body system [l| 
and offer unique opportunities to explore and under- 
stand its out-of-equilibrium behavior. Particularly inter- 
esting are one-dimensional systems which can be more 
or less closely approximated by integrable models and, 
therefore, allow one to study the effects of integrability 
and deviations from it. One-dimensional bosons [2| and 
fermions [3|, |j| are, in fact, quasi- ID systems in which 
the transverse motion is frozen out by a very tight con- 
finement. By integrating out the transverse motion in 
a two-body problem one obtains a ID delta-function po- 
tential for the effective interatomic interaction Q . In this 
manner quasi-lD spinless bosons and spin-1/2 fermions 
are modelled by the integrable Lieb-Liniger @ and Yang- 
Gaudin |7|, |8| hamiltonians respectively. 

Quantum integrability implies that the A^-body scat- 
tering is nondiffractive, i.e. the outgoing momenta of 
particles are restricted to be only rearrangements of the 
incoming ones [9|. The two-body scattering in ID is nec- 
essarily diffractionless for equal mass particles. However, 
already for N = 3 the nondiffractive rule has nontriv- 
ial consequences. In particular, the probability of re- 
flection and break-up in atom-dimer collisions vanishes 
independent of the collision energy [l^j • This manifestly 
quantum phenomenon, which has never been observed 
experimentally, is a microscopic analog of the reflection- 
less scattering of solitons. 

Small imperfections of a realistic quasi-lD system com- 
pared to its idealized integrable counterpart are practi- 
cally irrelevant as long as we are interested, for exam- 
ple, in the equation of state or thermodynamic quanti- 
ties. However, the statistics of energy levels, localization 
of eigenstates, decay of excitations, thermalization, re- 
sponse to an external perturbation, transport, and other 
dynamical properties are sensitive to deviations from in- 



tegrability. Muryshev et al. [Tl[ have found that virtual 
excitations of transverse modes are responsible for the 
dissipative dynamics of dark solitons in a weakly inter- 
acting quasi- ID Bose gas. They have shown (see also 
\l% [l3| ) that the effect of the transverse degrees of free- 
dom can be accounted for by adding a local three-body 
interaction term into the nonlinear Schrodinger equa- 
tion (equivalent in this case to the Lieb-Liniger model). 
Yurovsky et al. [ij] have numerically solved a purely 
ID three-boson problem near a Feshbaeh resonance and 
found that the integrability is lifted when the two-body 
coupling constant becomes energy-dependent. 

In this paper we solve the quasi- ID atom-dimer scat- 
teri ng p roblem by extending the approach of Mora et 
al. [15|, [la] to the case of finite collision momenta. In 
the nearly-integrable limit, \a\/lQ <^ 1, we analytically 
calculate the scattering amplitudes and find that at a 
finite collision momentum q the leading contributions 
to the reflection and break-up probabilities are propor- 
tional to aP' /l}^q^ for fermions and to a^ /l^q^ for bosons. 
We show that the diffraction in the fermionic case origi- 
nates mainly from the effective range corrections to the 
quasi-lD two-body scattering amplitude. In contrast, the 
bosonic result is due to a diffractive scattering of the 
atom by the transversely excited part of the dimer, con- 
sistent with the local three-body potential of Refs. [ill - 

The strong a-dependence of diffraction suggests a way 
of creating a system with fixed thermodynamic properties 
but with a tunable integrability parameter, the fermionic 
case being more practical due to the higher power in the 
a-dependence and, as we also discuss, due to the suppres- 
sion of inelastic processes. A four-fold increase of a and 
doubling of ^o does not change the effective ID model (the 
effective coupling constant is unchanged) , whereas the re- 
flection probability in atom-dimer collisions increases by 
a factor of 64 in the fermionic case and by a factor of 16 
for bosons. This effect can be verified by colliding clouds 



of atoms and dimcrs in a similar but more tunable fashion 
as the quantum Newton's cradle experiment of Kinoshita 
et al. [l4|- We point out that a quasi- ID two-component 
''"K mixture close to a zero crossing for the interspecies 
scattering length is a promising candidate for exploring 
this and other phenomena that are sensitive to deviations 
from integrability. 

The paper is organized as follows. In Sec. |ll] we in- 
troduce our notation and briefl y re derive the three-body 
integral equations of Refs. [1^, [l^. In Sec. HIT] we sepa- 
rate the "integrable" part from the "perturbation" and 
develop the corresponding perturbation theory by con- 
structing the resolvent of the integrable part, the small 
parameter being the ratio of the dimer binding energy to 
the confinement frequency. Then, we present the pertur- 
bative results for the reflection, transmission, and break- 
up probabilities as well as the three-body recombina- 
tion rate constant in the fermionic (Sec. IIV[) and bosonic 
(Sec.|V| cases. The relaxation and recombination to deep 
molecular states are discussed in Sec. |Vll In Sec. I VIII we 
propose an experiment in which one can control the de- 
gree of integrability and discuss effects that can be stud- 
ied there from a more general perspective. We conclude 
in Sec. IVini 



We now introduce an auxiliary function / proportional 
to the regular part of ^ in the vicinity of ri = 0: 

lim 9[ri*(ri,r2)]/9ri = -/(r2)/47ra. (6) 

ri— >0 

The function / can be considered as the wavefunction 
for the atom-dimcr relative motion. In particular, the 
atom-dimer scattering phase shifts can be read off of its 
long-distance asymptote. In our problem the atom and 
the dimer are initially in the transverse ground state and, 
since the total radial angular momentum is conserved, it 
is sufficient to consider f{x,p) = f{x,p). In the rest of 
this section we will derive an equation for this function 
and discuss its relation to the scattering amplitudes. 

Comparing Eqs. ([3]) and dH) we rewrite Eq. ([2|) in the 
form 

{Ho - -B)*(ri,r2) = ./(r2)<5(ri) - /(f2)<5(fi). (7) 

Then let us expand / in eigenstates of a single-particle 
quasi-lD Hamiltonian: 






(8) 



II. THE THREE-BODY INTEGRAL EQUATION 

Let us first consider the fermionic t-f-i system in a 
cylindrical harmonic trap. We use units h = m = lo = 1, 
where m is the particle mass and u is the frequency of the 
radial confinement. The oscillator length l^ = yjh/muj = 
1 is our unit of length. After separating the center of 
mass motion the noninteracting three-body Hamiltonian 
reads: 



Ho 



-p?/4 + P2/4-2, 



(1) 



where ri = {xi^pi] is the distance from one of the f- 
atoms to the |-atom, and -\/3r2/2 = {\/3.T2/2, \/3p2/2} 
is the distance from their center of mass to the second 
f-atom. The last term in Eq. ([T]) shifts the ground state 
energy to zero. Including interactions the Schrodingcr 
equation reads: 

(iJo--B)*(ri,r2) = -[t/(ri)+[/(|ri-V3r2|/2)]*(ri,r2), 

(2) 
where the interspecies interaction is taken as the zero- 
range Fermi pscudopotential 



U{r)- = ^TTa5{v)d{r-)/dr. 



(3) 



The wavefunction of the system should be antisym- 
metric with respect to the permutation of identical t- 
fermions, i.e. 



where 



*(ri,r2) = -*(ri ^fi,r2^f2), (4) 

fi = fi(ri,r2) = (ri - V3r2)/2, 

f2 = f2(ri,r2) = -(V3ri+r2)/2. (5) 



where Rm{p) = im(p^/2) exp(— p^/4)/-\/27r are normal- 
ized radially symmetric eigenfunctions of a 2D harmonic 
oscillator, 

(-V^ + pV4-l)i?m(/o)-2TOi?™(p), (9) 

and Lm are Laguerre polynomials. The solution of 
Eq. ([7]) can now be written as 



m=0 



dk 



vl'(ri,r2) = J2 / ^/™(fc) [^™(P2)e"=^=GB_ 



fc2-27 



Rm{P2)e^ ^'^GE-k^~27 



^{h)] 



.(ri) 
(10) 



The Green function Ge{j^) satisfies the equation (— V^ + 
p^/4:— 1 — E)Ge{tc) = S{r) and can be explicitly written 
as 

^^^") = Jo (4^)3/2 V7sinhr ^" 



1 

Airr 



C{l/2,-E/2) 
A-Ks/2 



o(r), 



(11) 



where we also present the two leading terms in its small-r 
expansion. In Eq. (Ilip C, is the Hurwitz zeta function. 

We now substitute pH]) into the left hand side 
of Eq. ([6]), multiply the resulting equation by 
-Rn(p2)exp(— ipa;2), and integrate it over v-i. The result 



C(l/2,pV2-£;/2 + n) 
47rV2 



Uv)' 



oo 

-E 



Mnrnfniip) 



Up) 

Ana 
(12) 



where the operator Mnm is defined as 



III. PERTURBATIVE FORMALISM 



Mnmfrnip) 



with the kernel 



rlk 

—M^„AE - (4/3)(p2 + k^ ^ pk)]f„,{k) 

(13) 



OO /"OO 



L„(z)L„(z/4) 



X exp 



/o ^0 
5 + 3 coth T 



cxp[(g-2m)r] 
'^^ l_exp(-2r) "^ ^ ^ 



The function /g (p) can be thought of as the atom-dimer 
scattering wavefunction in momentum space and its sin- 
gularities can be made explicit: 



Io{p) = 27r(5(p - q) - igo{p)2q/{p'^ 



iO), (15) 



where the scattering momentum q > satisfies 

C(l/2, <7V2 - E/2) = -V2/a = C(l/2, Ed/2). (16) 

Here we introduce the dimer binding energy Ed > 
which equals q^ ~ E and, therefore, the second equal- 
ity in Eq. ((T5)) directly relates a and Ed- In the case of 
small and negative a the binding energy is also small. 
Expanding the Hurwitz zeta function [see Eq. ((20)) ] we 
obtain S^ « a^ ^ 1. 

The function go{p) is smooth, and its values sd, p ^ q 
and p = —q are related to the atom-dimer transmission 
and reflection amplitudes; 



^l?) = 1+50(9), r{q) ^ go{-q). 



(17) 



The functions /„ for n > also have poles (on the real 
axis) if the collision energy is high enough to excite the 
relative atom-dimer motion to the n-th transverse state, 
i.e. q^ > 2n. The correct rule of integrating the poles 



can be enforced by the ansatz /„ (p) 
q'^ - iO), where g„ = y^ 



'2ig„ip)qn/{p 



2n and the functions gn{p) 
are smooth. The amplitudes of the forward and backward 
propagating waves in the n-th channel are given by gn (qn) 
and gni—Qn) respectively. 

The break-up channel opens for q^ > Ed, i.e. for posi- 
tive E. We can still use Eq. (|T0| in this case by choosing 
the retarded Green function Ge>o which requires that 
the flux of the released atoms be directed to infinity in 
the plane {xi,X2}. The retarded Green function is ob- 
tained by the analytic continuation of Ge<o to positive 
energies along a contour in the upper halfplane, or, equiv- 
alcntly, by substituting E -^ E + iO. One should then 
proceed with solving Eq. (|12p respecting the branch cuts 
and the rules of residue integrations. For example, \J —E 
should be substituted by —i\fE for £■ > 0. The break-up 
probability can then be calculated, for example, from the 
equation 

[9V2] 
n(g) = l-|t(<z)p-|r(g)p- Y. l5n(9n)P + l5n(-9n)P, 



n=\ 



(18) 



which follows from the atom number conservation law. 



Mora et al. ^M,\lM have derived Eq. p2|) . projected 
it to the lowest transverse channel by setting f„>o = 0, 
and solved it numerically at zero collision energy. In this 
manner they determine the ID scattering lengths for the 
even and odd channels, a^d and bad, as functions of Ed- 
They also note that in the limit Ed ^ Eq. (|T2|) takes 
the form of a purely ID integral equation, the solution of 
which can be found analytically from the Bethe ansatz. 
The transmission and reflection amplitudes are known in 
this case for any collision energy: 



i(0)(9) = 



-y/q^-E + VSiq 



V? 



E + V^iq 



A'>\q)=0, (19) 



and the break-up reaction probability also strictly van- 
ishes. 

In order to derive the corrections to the amplitudes 
([T9| we return to Eq. ([T2l) and separate out its integrable 
part in the following manner. We represent the Hurwitz 
zeta function as 



^1^'^) = V^ 



2 exp(£;r) 
exp(2T) 



1 



1 



dr 



where we make explicit the part diverging at £' — )• 0. 
Similarly, Eq. p^ for the ground transverse channels 
reads 



Moo{E) 



1 



7rV3 



1 



exp(i5r)dr 
4cxp(2t) 



1 



(21) 



Now, by using Eqs. (|T6|) . (|20)) . and (|2T|) we rewrite the 
n = part of Eq. ((HI) as 



{L - X,)Mp) - V{p)fo{p) + m'Mp) 



where the operator L is defined as 



Y Momfmip), 
(22) 



m>0 



Lfip) = 



fip) 



VS f dk 



Any/p^ - E 47r 



fik) 



27rfc2 



_ pz 



kp-SE/A' 
(23) 



\q = l/(47r-\/g2 — E), the function V{p) equals 

exp{ET — q^r) — exp{ET — p^r) dr 



1 



^(^) = -. 



exp(2T) 

and M' is an operator defined by 
1 



1 



ITT 

(24) 



M'fip) = 



VStt J 2tt 



?/(fc) 



dr- 



ET-{4/3){p^+k^+pk)7 



4e 



2t 



1 



(25) 



In the limit Ed ^ the right hand side of Eq. ([22)) 
vanishes. Neglecting it one arrives at an integral equation 
describing the purely one-dimensional integrable case |16l . 
Il8j which is our zeroth order starting point. Since the 



atom-dimer scattering solution is known from the Bethe 
ansatz for any momentum q G (—00,00), the operator L 
is easily diagonalized. Namely, the eigenstate equation 



LXijip) = hXqip) 



(26) 



is solved by 



Xq{p) = 27r5(p-g)-i 



t(0)(q)-l i[t(0)(g)-l](p + 2g) 



p — q — iO p^ + (f' + pq — 'SE/A 

(27) 
Accordingly, as the zcroth order solution of Eq. (|22|) we 
take 



f(0) 



Aiji Jm>0 



0. 



(28) 



Consider for simplicity that the collision energy is 
smaller or of order Ed-, so that there is a single small pa- 
rameter Ed <^ I. We then substitute ((281) into Eq. (p2)) 
and estimate the magnitude of different terms by using 
the fact that typical momenta involved in Xq{p) ^^^ of 
order y/Ed- We see that L ^ Xq ^ XjsfEd. Expanding 
Eq. ([M)) to the leading order in E, q, and p we get 



y(p)«[C(3/2)/16V27r](/ 



^\ 



(29) 



where C(3/2) sa 2.61 is the Riemann zeta function. From 
Eq. (|29| we see that V ^ Ed when acting on Xq{p) a-nd 
the leading order correction to /o, which can be written 

as {L - Xq)~^Vip)xq{p), is of order E^^^Xq- 

To understand the physics behind the F-term in 
Eq. (|22|) imagine a purely- ID problem in which f and 
4, fcrmions interact via a delta-function potential, the 
strength of which depends on their collision energy. If 
this dependence is chosen such that the corresponding 
scattering amplitude matches the one obtained for two 
quasi-lD atoms [1^, we arrive exactly at the three-body 
equation {L — Xq)f{p) = V{p)f{p). The V-ieruv then 
reflects the effective range corrections to the two-body 
interaction due to virtual transverse excitations. That 
this type of perturbation breaks integrability in a system 
of three identical ID bosons has been shown numerically 
by Yurovsky et al. |14| . 

As a side remark we note that so far we have been 
considering the idealized zero-range 3D interaction po- 
tential ([3]). Corrections corresponding to the finite range 
of a realistic 3D potential or to a finite Feshbach res- 
onance width can be incorporated into the formalism of 
Sec.|ll]by introducing an energy dependent 3D scattering 
length, 1/a — 7- l/a(tcoii)- The collision energy ecoii is de- 
fined as the kinetic energy, — V^^, of the relative motion 
of two atoms when they are close to each other but still 
outside of the support of the 3D interaction potential. If 
their relative motion with respect to the third atom is 
described by the wavefunction Rm{p2)e^^^'^ one can see 
from Eq. dJ) that ecoii ^ E - k"^ + I - 2m. Adding the 
effective range term, 1/a — > 1/a— (ro/2)ecoH, one arrives 
at Eq. (dl]) where 



(30) 



The right hand side of Eq. ((22)) then acquires an addi- 
tional term 



V{p)Mp) = -(ro/87r)(/ - q')fo{p), 



(31) 



which has the same form as Eq. ([29)) . Therefore, if we 
know the correction to / due to the perturbation V, the 
perturbation V does not require any special treatment. 
We simply multiply the F-result by 1 — 2-\/2?'o/C(3/2). 
Note that V '^ V only when rg is comparable to the 
transverse oscillator length, which is typically of order 
100 nm and much larger than the physical range of van 
der Waals potentials for neutral atoms. Therefore, one 
should care about V only when tq is anomalously large, 
in particular, in the case of a very narrow Feshbach res- 
onance. 

Let us now turn to the second term in the 
right hand side of Eq. (|^^ . Approximating 

^Er-{4/3)ip^+k'+pk)r ^^ -^ -^^ gq_ ^ ^^ obtain M'fip) w 

log(2/-\/3)/\/37r/ /(p)(ip/27r, i.e. this term is local in 
the position representation. Adopting the confinement- 
induced resonance terminology the M'-term describes 
the interaction of the third atom with the "closed chan- 
nel" molecule formed by the first two atoms. If we were 
dealing with bosons (see Sec. |V| the M'-term would give 
a correction to /o of order EdXq , i-e. it would be more 
important than the V-term. However, for the fermionic 
f-f-i system J Xq{p)dp = which follows from the fact 
that three atoms can not be at one point in space. In 
this case the leading contribution to the M'-term is ob- 
tained by further expanding Eq. (|25)) and corresponds 
to the odd-channel interaction of the "closed channel" 
molecule with an atom. We will quantify this interac- 
tion later [see Eq. (|57)) ]. Now it is sufficient for us to say 

that Af ~ eJ when acting on Xqip) and is thus less 
important than the y-term. 

Finally, in order to estimate the contribution of the 
last term in Eq. (P^ we have to consider Eq. (|T^ for 
n > 0. To the leading order we obtain 



Up) « -{l/Xq)MnOXq{p)- 



(32) 



1/a -^ 1/a - {ro/2){E - p^ + 1 - 2n). 



From Eq. ([13]) it follows that M„o(^ ^ 0) = const and 
similarly to the operator M' the order of magnitude of 
/„ depends on whether J Xqip)dp vanishes or not. For 
fermions it does and we have /„ ^ ^'dXq- Substituting 
this result into the last term in Eq. ()22p we see that 
the contribution of the higher transverse channels can be 
safely neglected even compared to the M'-term (cf. [16|). 
This is what we will do and until further notice we omit 
the subscript of the function fo{p). 

In the first approximation Eq. ([^^ can be solved by 
substituting f^'^' in its right hand side and by inverting 
the operator L — Xq. We keep the operator M' in play 
because its leading order contribution is still larger than 
higher order terms related to V. Accordingly, we call 
/(I) = (L- Xq)-^iV + M')f^°'> the first order correc- 
tion to / although it actually contains the leading or- 



der term proportional to eJ Xq a-s well as the next-to- 
leading term c>c E^Xg- To diagonalize the operator L 
we expand /(p) in the basis of its eigenstates Xg [see 
Eq. ([2B1)]: 



/(P) 



akXk{p)dk/2TT. 



(33) 



One can directly show that these eigenfunctions are or- 
thonormal in the sense 



(91 1 92) 



Xqi (P)X92 {p)dp/2'K = 27r5{qi - 92)- (34) 



Here Xg = Xq if ^ < 9^: which is always the case be- 
low the break-up threshold. Otherwise, as mentioned in 
the end of Sec. HI] we require all the functions that we 
are dealing with be analytic in the upper halfplane of 
complex variable E. In particular, the functions Xg ^nd 
Xg are obtained by analytic continuation of (P7)) and its 
complex conjugate from the E < q^ to E > q^ part of 
the real axis following a path in the upper halfplane. 

Changing the basis from plane waves to Xg diagonalizes 
the left hand side of Eq. 



(Afc - Xg)ak 



\V + M'\k')ak'dk'/2TT. 



(35) 



The matrix elements in Eq. psp are 



{k\V\k') 



C {k\/k'^ - E - fcVfc2 - E)/{k - k') 



in 



4(fc2 + /s'2 + kk')/3 - E 



X {^/k^ - E + iVSk){\/k'^ - E- iVSk') (36) 



and 



C" 



{k\M'\k') w —{iVSk + Vk^-E){iV3k' - Vk'^ - E), 
in 

(37) 
where 

C = \/3C(3/2)/8\/2 « 0.400, 
C" = Li2(l/4)/\/3w 0.155, 

and Li2(l/4) is the polylogarithm function. When calcu- 
lating the matrix element p6p we use the approximation 
(P^ and in Eq. ([57)) we also retain only the first nonvan- 
ishing term in the expansion of Eq. p5|) . Equations (|35f 
[37| are therefore vahd for q'^,k'^,k'^,E <t: 1. 

The iterative solution of Eq. ([55]) is now straight- 
forward. We substitute the zeroth order term a 
2nS{k — q) into its right hand side and obtain 



(0) 



a['^^{k\V + M'\q)/i\k~X,). 



(38) 



We then find f''^\p) from Eq. ((33)) . The corresponding 
corrections to the refiection and transmission amplitudes 
(fro)) can be calculated from the residues of f'^^\p) at 
p = —q and p = q, respectively [see Eqs. (|T5)) and (flT)) ]. 



These residues are obtained by performing integration 
in Eq. (|55| in the vicinities of fc = =Fg. Close to these 
points the matrix elements ([36]) and (|37| are smooth and 
l/(Afc - Xq) « -87r(g2 _ Ef/^/{k^ - q^ - iO). Here the 
positions of the poles with respect to the real axis are cho- 
sen such that there is no incoming wave with momentum 
—q. The corresponding correction to the transmission 
amplitude reads 



t^^\q)/t^°'>{q) = ~4ni{q^ - Ef/^{q\V + M'\q)/q 
.^^ C{l-q'/Ea)-C'El'\l + ?,q^/Ea) 

and for the reflection amplitude we get 

r(i)(9) = -Ani{q^ - Ef'^{-q\V + M'\q)/q 



(39) 



1/2 



i ^^''^"/"^-'-^'^'- (i-../v7i;)',(4o) 



-lE: 



Equations p9p and (|40)) give the first two leading correc- 
tions (ex E"^ and oc EJ ) to the atom-dimer transmission 
and reflection amplitudes for small E^. The validity of 
these equations requires that the scattering momentum 
be in the interval i?J ^ g <C 1. In particular, the atom- 
dimer collision energy q^ can be above or below the break- 
up threshold, Ed- The unphysical divergence of t'^^\q) 
and r'^^-*(g) at very small momenta q "^ E^ <^ \fEd is 
a consequence of the fact that in ID any weak interac- 
tion becomes strong at sufficiently low energies and the 
Born approximation (which is our first iteration) neces- 
sarily fails. To illustrate this we note that Eq. ((35)) at 
low energies (fc, k' ,q -^ y/Ed) reads 



where 



(fc2 _ q^)ak + 2Q 



,3/2 



ak'dk'/2n:^0, 



(41) 



= 4nE^J'{0\V + M'|0) = EJ{C - C yj Ed) (42) 



is a small characteristic momentum. Equation (|4ip is 
nothing else than the ID Schrodinger equation in momen- 
tum space describing the scattering on a weak (5-function 
potential. Its exact solution is 



Ctk 



2nS{k ~q)- 



2q 



g + iQ fc2 — g2 



tO 



(43) 



We sec that a^ is appreciable only at small momenta 
^ q <^ \/Ed, which justifies Eq. (01]). Note that solv- 
ing Eq. ((4ip in the first Born approximation we would 
miss the term iQ in the denominator of Eq. (H5)) . which 
would lead to the 1/q- divergence, exactly as we observe 
in Eqs. ^ and (gD)). 

Now substituting Eq. ((43)) into Eq. ((33)) and calculating 
the residues of f{p) at p = ±q we obtain the transmission 
and reflection amplitudes: 

r{q) = t{q)/t^'>Hq) - 1 - -lQ/{q + zQ). (44) 



Equations (|4m44p are valid for q^ <^ Ed- This range of 
collision energies has a large overlap with the interval of 
validity of Eqs. (p9| and (|40| . and they can be easily 
matched with the low-energy result (|44p . A universal 
result can be obtained simply by making the substitution 
q^f q + iQ m the denominators of Eqs. (f39| and pO]). 

Equation (^^ demonstrates that the integrable limit 
-Ed — > (Q — > 0) and the limit g — >■ do not commute. 
In the near-integrable case {Ed <C 1) the reflection prob- 
ability is of order 1 for collision energies 9^ ^ Q^ oc _Ej. 
Outside of this small region |r((7)p ex E\lq^ . Recalling 
that in the near-integrable case Ed ~ a^, we observe a 
very strong decrease of the reflection probability with de- 
creasing \a\. Namely, at a fixed collision energy |rp ex a®. 

The even, Fs{q), and odd, Fp{q), scattering amplitudes 
are related to the transmission and reflection amplitudes 
by the equations 

FM = [t{q) + r{q) - l]/2, Fp(g) = [t{q) - r{q) - l]/2, 

(45) 
and thus can be calculated to the first order from 
Eqs. ([M)) and PO]) for any collision energy. The even, aad, 
and odd, bad, scattering lengths can be defined through 
the effective range expansions of Fs and Fp valid for 
q^ « Ed-. 



Fs{q) 



1 + 2iaadq/Vi + 



and 



Fp{q) 



V^z/(26adg) + l-2zep9/\/3- 



(46) 



(47) 



where the rescaling of q is due to the fact that the atom- 
dimer scattering momentum in our units equals 2(7/a/3, 
which is the Fourier conjugate of the atom-dimer distance 
^/ix/2. In Eq. (|47)) we also introduce the odd-channel 
effective range ^p, so that the expansions orders of Fp 
and Fs match. The even scattering length equals (cf. 

M) 

aad = ^E^ ^^Ed + oiE^' ), (48) 

the odd scattering length and effective range are given 
by 

-L=^-^EJ- %eT + o{eT) - 4^' (49) 

Kd V3 V3 Vs 



and 



\e-^" - "^Ed + 3V3C'Er + oiET). 



(50) 

Let us now say a few words about higher order terms. 
The second order correction to a^ reads 



a 



(2) 



1 



Afc — A„ 



dk' {k\V + M'\k'){k'\V + M'\q) 



27r 



At' — A„ 



The integral in Eq. (|5T|) diverges at large k' if one uses 
the formulas p6ll37p . Investigating Eqs. (|24j|25p we see 
that the exact matrix elements start decaying much more 
rapidly than the approximate ones p6ll37p at momenta > 

I. Therefore, in order to estimate a\, we can introduce 
a cut-off at fc' -- 1 in Eq. ^^ and still use Eqs. (|36ll37p . 
One can then directly show that the leading term in a), 
is by a factor ~ ^/Ed smaller than the M '-contribution 
to a^, . This justifies our keeping the operator M' when 
calculating the first order terms. 



IV. REACTION PROBABILITIES 

We can now discuss the reflection, transmission, and 
break-up probabilities, |?'((7)p, |i(q)p, and F;,(g). The re- 
flection probability is approximated by \A^' (9)^ ensuring 
the two leading terms: 



\r{q)?^C 



.,EJ f Ed + 3q^ 
q^ \Ed + qy3 



1-2 



C 3Ed + q^ 

C i^/Wd 



(52) 

The break-up probability can be obtained from the equa- 
tion Pb(q) = I - |r((7)|2 - |t((7)p, in which \t{q)\'^ is cal- 
culated to the same order as |r((7)p in Eq. ([52]): 



n(.)^H^^^H.)PHt^^H.)P^2Re ^"^(^J+;p^) . 

(53) 
Here v'^'{q) is the correction to the transmission ampli- 
tude derived from Eq. (|5ip in the same manner as t^^'{q) 
is derived from Eq. ([38]): 



(51) 



^(^)_ 47rz(g^-£:)3/2 r dk' {q\V + M'\k'){k'\V + M'\q) 

T^" q J 2^ Ag - Afc' ■ 

(54) 
From Eq. ^ one can see that Re[t'^^\q)/t^°^q)] = 0. In 
fact, if Pbiq) contained first order terms, we could play 
with the sign of the corresponding perturbation operator 
and make the probability negative, which is not possi- 
ble in principle. Therefore, the operators that we have 
neglected earlier in favor of T^^ + M' would contribute to 
Eq. ([53]) only to the second order and thus we do not 
exceed accuracy by keeping the second order term t'^^ 
Anyway, the need for the second order terms here is for- 
mal as we know that the break-up probability can in 
principle be derived from the first order solution af, by 
restoring the wavefunction ^(ri,r2) and calculating the 
outgoing flux of free atoms. 

One can see that the real part of t^^V^^^'* in Eq. ([5i)) 
originates only from the integration interval k' S 
{—^/E, y/E) and from the residues at k' = ±q. The con- 
tribution of the latter cancels the first two terms in the 
right hand side of Eq. ([55]) . This means that Pb{q) = 
below the break-up threshold as it should. For E > we 



have explicitly 

Pb{q) -= — / Im 

q J-ve 



{q\V + M'\k'){k'\V + Ar\ql^^, 



Xk' — Xq 



and performing the integration we arrive at 



(55) 



Pbiq) 



C^El{q^-Eaf 
q^iEd + qysr 



^C El 



IQEdq^li + q^ 



C {^/E2+\q\f 



(56) 

Just above the break-up threshold Pi, oc E^. This is due 
to the ferniionization of unbound states at low energies - 
the probability to find three unbound interacting atoms 
close to each other scales with E^ as in the case of three 
identical fermions. 

Quantitative comparison of Eqs. ([5^ and ([55]) shows 
that |rp is significantly larger than Pf, even well above 
the break-up threshold (in the extreme limit Ed ^ the 
two curves intersect at q^ /Ed ~ 12), and, for example, 
for the collision energy q^ = 2Ed the ratio \r\^ /Pt ^ 40 
for Ed ^ 0.1. This means that the dissipative dynam- 
ics of a dimer immersed in a gas of atoms is dominated 
by the reflection rather than the break-up. On the other 
hand, the latter is a chemical reaction and its rate can be 
measured by monitoring the evolution of the population 
of dimers. The same holds for the break-up's inverse, 
the three-body recombination. The three-body recom- 
bination rate determines, for instance, the rate of dimer 
formation in a super- Tonks state - the state of the sys- 
tem with negative coupling constant but without dimers 
[20l - [2a | (the bosonic super- Tonks gas has been recently 
observed [H]). 

The three-body recombination rate constant a-f^i(E) 
is readily obtained from Pb{q) by using the principle of 
detailed balance. Indeed, the three-body recombination 
rate for a single f-f-i triple per unit length in the center 
of mass reference frame equals 2a^^^{E). Here the factor 
2 comes from the fact that the product of densities riin:^ 
is twice the number of t"t"4- triples per unit length. This 
rate multiplied by the density of unbound states (of a sin- 
gle triple) p^^i{E) = l/(4-\/37r) should equal the break-up 
rate Pb('z)V3'Z multiplied by the density of atom-dimer 
states Pad{E) ~ l/(-\/37rq), where E = —Ed + q^ ■ In this 
manner we get 

«tn(^) = «m(^) - 2\/3Pb(g), (57) 

which should be multiplied by h/m in order to restore 
the dimensions. We should mention that the rate con- 
stant (|57p is averaged over a uniform distribution of the 
initial three-body unbound states in a small energy in- 
terval close to E, i.e. in the ergodic approximation. For 
highly nonthermal distributions one should speak about 
the differential recombination rate of a particular ini- 
tial unbound state characterized by the asymptotic mo- 
menta (rapidities) and an auxiliary index related to the 
parametrization of the f-f-j wavefunction in the nested 
Bethe ansatz picture [2J]. For each such state one can 



determine the zeroth order function f^^\ expand it in 
the basis of Xk, substitute the corresponding coefficients 



.,(") 



into the right hand side of Eq. ([55)1 , determine a 



(1) 

k ' 



and finally derive the outgoing atom-dimer flux from the 
residues of the poles of f^^\p) at p = iq. However, 
Eq. ([57)1 holds for a thermal gas (in this case E ^ T) and 
is also useful for systems in which the atomic momentum 
distribution is not extremely exotic. In particular, we be- 
lieve that Eq. ([57]) with E ^ Ep gives a good estimate of 
the recombination rate constant in a degenerate super- 
Tonks gas, provided the Fermi energy Ep = -n^r? 12 is 
much smaller than the dimer binding energy Ed- 

Remarkable is that in the case E ^ Ed the recombi- 
nation rate constant (j57p is independent of Ed^ i.e. in- 
dependent of the interatomic interaction strength. The 
reason for this is the following. On the one hand we have 
a nonintegrable perturbation, which acts on three atoms 
when they are close to each other. The squared modulus 
of this perturbation (relative to the zeroth order terms) 
scales as E\ and it does vanish in the integrable case. On 
the other hand the local three-body correlation function 
(probability to find three atoms close together) is propor- 
tional to {E/EdY- We see that when we multiply these 
two factors the i^Jfj-depcndence drops out. In a thermal 
gas we have a^^^ ex {h/rn){T/huj)^. 



V. BOSONIC CASE 

In the case of three identical quasi-lD bosons the 
derivations of Sees. |TT] and IIIII are essentially the same. 
The modifications are related only to the facts that the 
bosonic wavefunction is symmetric and that all three 
atoms interact with each other. Accordingly, the bosonic 
version of Eq. (|12l) differs from the fermionic one by an 
extra factor —2 in front of the operators M„m [lal- The 
operators M' and Afom in Eq. (|22p and the integral in 
Eq. (P5|) should also be multiplied by this factor. This 
changes the properties of the eigenstates and the struc- 
ture of the spectrum of the operator L (we now have a 
trimer state separated from the continuum). 

In the bosonic case J Xqip)dp is finite since three 
bosons can be at one point in space. This fact leads to im- 
portant qualitative differences in between the fermionic 
and bosonic cases. Analyzing different terms in the 
bosonic version of Eq. (|22p in the same manner as we 
did for fermions in Sec. IIIII we see that the M'-term is 
the leading perturbation: 



M'fip) 



log(4/3) CM 
2V37r J 27r 



/(fc). 



(58) 



It is of order M' ^ yEd when acting on the bosonic 
Xq{p)j whereas V is still of order Ed as in the fermionic 
case. One can also easily show that the last term in 
Eq. (1221) which represents the coupling to higher trans- 
verse states is of the same order of magnitude as the V- 
term. Therefore, the next to leading order correction is 



more difficult to obtain compared to the fermionic case. 
Here we restrict ourselves to the leading order correction 
originating from the M'-term, but before presenting the 
results let us briefly mention the effect of a narrow Fes- 
hbach resonance (cf. [1J|)- In this case we should add 
the term V given by Eq. (|3T|) to the right hand side of 
Eq. ([22|) . We see, however, that it is larger than the 
A'/'-term only when the effective range tq is larger than 
the size of the ID dimer l/^/E^, which is an extremely 
restrictive condition. 

Formally, the il/'-term, local in real space, can be un- 
derstood as a modification of the one-dimensional inter- 
action between two bosons in the presence of a nearby 
third boson. One can show that calculating the correc- 
tion associated to this perturbation is equivalent to solv- 
ing the Schrodinger equation 

{-^U,^2} - 2VE~d[S{x,) + E± S{x,/2 ± V3X-2/2)] 
-8V3\og{A/3)EdS{xi)5{x2) - E^ *(xi, xa) = (59) 



to the first order in the perturbation 
-8V3\og{A/3)EdS{xi)d{x2) (cf. [iMl)- The un- 
perturbed operator in the first line of Eq. (j59|) is 
diagonalized by using the Bethe ansatz. The standard 
first order perturbation theory then gives the following 
results. 

The first order correction to the boson-diboson trans- 
mission amplitude reads 



t^'\q) 



A\ogi4/3)El/^l + 3qyEd 



t^"\q), (60) 



V3 q l + qySE/ 

where the zeroth order transmission amplitude equals 



f + VSiq/VE^ 1 + ig/ys^' 

The zeroth order reflection amplitude vanishes and in the 
first order we get 



.(1) 



(9) 



. 41og(4/3)g,^/^ 
V3 q 



' 1 - Vitq/VE^^ 



(62) 



Accordingly, to the leading order the rcfiection proba- 
bility equals Ir^^^^p. The break-up probability above the 
break-up threshold {E = q^ — Ed > 0) is given by 



161og^(4/3)£;^'/' jq^-Ed)^ l + Sq^Ed 
3 g2 (1^1 



^'^'^^•- - "2 "- ■ /Wdfl + qyiEd 

(63) 
and the three-body recombination rate constant equals 
arec{E) = 2-\/3Pb(q) as in the fermionic case. Here we 
also observe the scaling P^ oc arec oc E^ just above the 
threshold. This fact is related to the suppression of the 
local three-body density-density correlation function at 
low energies in an interacting ID gas. That the two- 
body correlator is also suppressed can be seen by looking 



at the differential rate. Indeed, the recombination rate 
for three bosons on the length L in a state parameterized 
by the set of rapidities {fci, fe, fca} equals 

J^rec(fcl,fc2,fc3) « mVllog'' {A/ i)L~^ [E'' / Ed) sixv^ {icj)) , 

(64) 
where = arctan[v3(fc2 — fc3)/(2fci — ^2 — k-^)] and we 
assume that the energy in the center of mass reference 
frame E = J^i^f/'^ ^ (X]i^"i)^/6 is much smaller than 
Ed- We observe that Vrec oc {ki — kj)^ for any pair of 
rapidities, if they are close to each other, {ki — kj)"^ <C E. 
Due to the fact that the nonintegrable perturbation 
scales with lower power of Ed compared to the fermionic 
case the refiection and break-up probabilities, as well 
as the rate of three-body recombination in the bosonic 
case contains one power of Ed less. In particular, the 
three-body recombination rate actually increases with 
decreasing Ed (although Ed should be kept larger than 
E). This counterintuitive phenomenon is explained by 
the fact that the local three-body correlator is propor- 
tional to (E/Ed)^, whereas the squared modulus of the 
perturbation is oc E"^. 



VI. RELAXATION TO DEEP STATES 

Let us compare the reflection probability to the prob- 
ability of relaxation to deeply bound molecular states. 
This inelastic process is extremely local. It takes place 
at distances ~ i?e ^ 1, where Re is the van der Waals 
range of the interatomic potential. The relaxation prob- 
ability is thus proportional to the probability of flnding 
three atoms at distances of order the oscillator length 
(unit in our case) multiplied by the recombination rate 
for three atoms conflned to a unit 3D volume. We stay in 
the near-integrable regime where |a| ^ 1 and Ed ~ a^. 
This gives the following results. 

In the fermionic case the probability of relaxation in 
an atom-dimer collision at the collision energy q^ equals 

Preliq) OC q-'E'J\l + 3qyEd){Re/\a\f+''', (65) 



where 7 « -0.2273 [25|. Comparing Eqs. (|65]) and ^ 
we see that the relaxation probability is always much 
smaller than the reflection probability. One can also show 
that the recombination to deep molecular states is much 
slower than the formation of shallow dimers in a gas of 
unbound fermionic atoms. Namely, for E < Ed the ratio 
of the corresponding recombination constants scales as 
ad,m/«m « EdiRe/\a\)^+^^ « 1. 

The relaxation probability in boson-diboson collisions 
can be calculated including the prcfactor. Indeed, in the 
zeroth order the wavefunction of an atom and a dimer 
normalized per (axial) length L is given by the axial 
Bethe ansatz wavefunction multiplied by the wavefunc- 
tions of the radial ground states. When the three atoms 
are at distances much smaller than 1 (the radial oscilla- 
tor length) this wavefunction is approximately constant 



and its modulus square equals |^P = {l/9Tr'^L)^/Ed{l + 
2>q^ /Ed)/{1 + q^/3Ed). Here we use the normalization 
corresponding to the "normal" laboratory coordinates 
in which the atom-dimer distance is not rescaled (we 
mean the factor V3/2 in the definition of r2). The re- 
laxation rate in this case equals Vrei = (Xdia < 0)|^p, 
where ad(a < 0) is the relaxation rate constant for 
cold thermal bosons in the uniform space. This quan- 
tity has been calculated numerically [2y, [23] and an- 
alytically [23 in the zero-range theory as a function 
of a, the three-body parameter, and the elasticity pa- 
rameter. We can write ad{a < 0) « Cda^, where 
Cd is a log-periodic function of a bounded from below: 
Cd > 1287r2(47r-3V3)coth(7rso)tanh7?* « lOV- Here 
So ~ 1.00624 and the elasticity parameter is typically 
77* > 0.1. Assuming the case Cd ~ 10'^ we have 



Preliq) = {L/V3q)Vrel 



6.5 '^ 



Sq^Ed 



q l + qySEd' 



(66) 



We observe that the ratio Prei/kP ^ ^5q/^/Ed, which 
means that the relaxation is the dominant process in 
boson-diboson collisions unless we find a Feshbach res- 
onance with a very low elasticity parameter and/or go to 
extremely low collision energies. 

The above analysis can also be performed for un- 
bound states of three quasi-lD bosons. Namely, the 
relaxation rate of a Bethe ansatz state with rapidities 
{fci,fc2,fc3} equals Vd{ki,k2,k^) = o-difl < 0)|*p, where 
|*P = {l/8n^L'^){E/Edfsm'^{3cj)). This rate scales as 
I /Ed ~ l/a^ in the same manner as the rate of three- 
body recombination to the weakly bound state (|M)) . In 
fact, these rates are related by Vd/vrec = 2.5 x IQ^^Cd 
and we see that the recombination to deep states actually 
dominates, assuming, for example, the value Cd ~ lO''. 



VII. DISCUSSION 

A direct consequence of the relaxation analysis of the 
previous section is that for bosons, at least in the near- 
integrable case (|a| <C 1), the thermalization due to the 
local three-body coupling of different Bethe ansatz states 
is likely to be slower than the relaxation and/or recombi- 
nation to deeply bound molecules. This means that such 
a gas is stuck in a nonthermal state during all its lifetime 
(which may be long because of the ID fermionization). 
In contrast, the relaxation processes in the fermionic case 
are suppressed and one should be able to observe the 
diffraction of momenta in atom-dimer collisions or the 
three-body recombination to a weakly bound state well 
before the gas decays. 

The dynamics of a mobile impurity in a ID gas is at 
present a very attractive theme of theoretical (s ee [291 ISOj 
and references therein) and experimental [3l[, |32[ stud- 
ies. The role of integrability in this problem has been 
discussed [33 but is still far from being well understood. 
The present paper prepares grounds for studies in this 



direction. A quasi-lD ti dimcr immersed in a Fermi sea 
of t species is a long-lived system in which, keeping the 
thermodynamic properties unchanged, one can modify 
the degree of integrability. Indeed, the thermodynamic 
quantities depend on the ID 'f-^ coupling constant, which 
is proportional to the product of the scatteringlength a 
and the frequency a; of the radial confinement [5| . Keep- 
ing this product constant we can still modify the atom- 
dimer reflection probability in a wide range by using the 
fact that its scaling with a and lo is drastically different. 
Namely, at a fixed collision energy we have |rp ex a^w^ 
[see Eq. dSH)]. 

Having established control over the degree of inte- 
grability one can study various physical problems, e.g. 
whether or not a dimcr looses its momentum while mov- 
ing through a gas of atoms, i.e. whether or not there is a 
friction force on the dimer. Based on our understanding 
of the atom-dimer scattering problem we can conjecture 
what happens if one collides a gas of atoms with a gas 
of dimers in a quasi-lD trap with a weak axial parabolic 
confinement. In the integrable limit these two clouds 
pass through each other without rcfiection and one ex- 
pects undamped relative oscillations. In contrast, as we 
increase the refiection probability these oscillations be- 
come damped and we expect to see none of them when 
|r((7)p > 1/A^, where N is the particle number. This sce- 
nario remains based on the microscopic few-body analy- 
sis and wc do not exclude nontrivial many-body effects, 
especially in the degenerate regime. Wc think that this 
problem deserves further experimental and theoretical in- 
vestigation. 

Note that the limit a — ?> does not at all correspond 
to the noninteracting case. Indeed, one can imagine 
a dimer in a state given by a wavepacket localized in 
momentum space around k and in real space around 
xq 4- fct/2. Let us assume that it passes an atom which 
is in a similar state but with fc = 0. This picture 
describes the relative motion of a dimer and an atom 
at high energies close to the bottom of the axial har- 
monic trap. After their collision the dimer wavepacket 
is centered around x^ + 5x + fci/2, i.e. its trajectory is 
shifted by 5x (the atomic wavepacket is then centered 
at —25x). By using the fermionic atom-dimer scatter- 
ing solution ([?f| and the transmission amplitude p^ 
with q = fc/2v3 one can show that this shift equals 
5x — y/Ed/{Ed + k"^ /A). For the boson-diboson scatter- 
ing we have fe = V;E;^[(£:d+fcV4)-i+(3£:d+A:Vl2)-i]. 
Note that if k"^ is of order Ed, the shift is proportional 
to Xj \fEd oc l/|a| and actually increases with decreasing 
\a\. The dimer passes the atom faster than it does in the 
noninteracting case (as if the atom-dimer interaction is 
attractive). If we now assume that there is a gas of N 
atoms with A^ ^ 1, their effect on the dimer's trajectory 
and on the frequency of its oscillations is appreciable. 
Note also that for a longitudinally trapped gas this shift 
is another source of nonintegrability even in the purely 
ID case [H. 

As a more concrete experimental proposal for investi- 
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gating the above phenomena we envision a setup, which 
is a combination of the fermionic experiment of Moritz 
et al. [3] on quasi-lD dimers and the bosonic quantum 
Newton's cradle experiment of Ref. [13| • Namely, one can 
follow the routine of the former and create a quasi-lD 
spin mixture (in our case spin-imbalanced) of two differ- 
ent hyperfine states of '^'^K close to a zero crossing for 
the interspecies scattering length and form molecules by 
adiabatically decreasing the scattering length from zero 
to a finite a < 0. One can then separate the dimers from 
atoms in momentum space by applying a Bragg pulse as 
has been demonstrated by Veeravalli et al. [35[ in the case 
of a 3D spin-mixture of ^Li. The transfer of ''^K dimers 
into a state with momentum 47r/A, where A = 767 nm, 
can be done by using counterpropagating Bragg beams 
with the frequency detuning i/2m = 16.97 kHz. Since the 
atoms remain at rest the atom-dimer scattering energy 
then equals i^2m/3 = 5.66 kHz. Assuming the radial con- 
finement frequency j/j^ = 100 kHz, which corresponds to 
the oscillator length ^o = 50 nm, this sets q^ = 5.66 x 10~^ 
in our formulas. For a = —25 nm the dimer size equals 
155 nm, its binding energy Ed « 10.5 x 10~^ (10.5 kHz 
in the laboratory units), and the reflection probability 
equals |r((7)p = 1.2 x 10~^. Under these conditions a 
cloud of, say, 30 dimers will pass through a cloud of 
the same number of atoms several tens of times with- 
out thermalization. Increasing a by a factor of 2 leads 
to |?'((7)p = 0.01 and we expect to see only very few os- 
cillations. Although the value of a in this case equals Iq, 
we show below that the binding energy calculated from 
Eq. P^ . Ed ~ 0.22, is sufficiently small to apply our 
perturbation theory. 
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FIG. 1: The quantity 1 — |t(g)p versus q'^ /Ed below the break- 
up threshold for several values of Ed- Numerical solutions of 
Eq. (|12|) are obtained by taking into account the first three 
(solid), the first two (dash-dotted), or only the ground (dot- 
ted) transverse states. The dashed lines is our perturbative 
result valid for small Ed- The dip at q^ /Ed = 0.5 in the case 
-Ed = 4 is discussed in the text. 



For this type of experiments the relevant values of the 
refiection probability are small, of order 1/iV, meaning 
that we can use the perturbation theory valid for Ed <C 



1. In order to see the limitations of the perturbative 
approach we have also numerically solved Eq. p^ and 
calculated |r((7)p and |i(g)P for various Ed- In Fig. [1] 
we show the quantity 1 — |t(g)P versus q^ /Ed below the 
break-up threshold for six different values of Ed from 
0.125 to 4. The perturbative result is shown as dashed 
lines and is given by Eq. ([52|) multiplied by q'^ / {q'^ + Q^) , 
where Q^ = C^EJ{1 - 2C'^/E^/C) [see discussion after 
Eq. (|44|) ]. We present it only for the three lowest values 
of Ed since for larger Ed the perturbation theory breaks 
down as expected - the next-to-leading order term in 
Eq. ([5^ becomes comparable to the leading one. 

For each value of Ed in Fig. [T] we also present three nu- 
merical curves: the solid lines are calculated by solving 
Eq. p^ projected to the first three transverse channels, 
i.e. we set fn>2 = 0, the dash-dotted lines correspond 
to the first two channels, /„>i = 0, and the dotted lines 
are obtained by projecting Eq. ([T2|) to the ground trans- 
verse state (/„>o = 0). For Ed < 0.5 the three lines 
practically coincide. For larger Ed the projection to the 
ground transverse state is insufficient, the dotted lines are 
far off. However, we clearly observe a fast convergence 
with increasing the number of kept transverse states. 

In the case Ed = 4: the point q^ /Ed = 0.5 (g^ = 2) 
is the threshold for the excitation of the relative atom- 
dimer motion to the first transversely excited state (see 
the end of Scc.[TT]). The corresponding branch-cut singu- 
larities are visible only in the solid and dash-dotted lines 
since the dotted line ignores all transverse channels other 
than the ground. In Fig. [T] the difference in between the 
quantities |r((7)p and 1 — |i((z)P exists only in the case 
Ed = 4: for q'^/Ed > 0.5 where there is a finite probabil- 
ity for the atom-dimer pair to be excited into the first 
transverse channel. 

The small-g expansion of the transmission probability 
starts with the term |i((7)P ~ {4:/3){aad + bad)^q^ which 
can be derived from Eqs. (|45II47[) . In fact, one can show 
that the coefficient in front of q"* is also proportional to 
o-ad + bad- Thcrcforc, the point at which aad — —bad is 
rather peculiar. In this case there is a large region of 
collision energies where atoms and dimers can be consid- 
ered impenetrable. Then a dimer immersed in a gas of 
atoms can move only by shoving the atoms on its way, 
which leads to its diverging effective mass, subdiffusive 
propagation dynamics, etc. (see, for example, [2^). This 
regime is opposite to the integrablc limit where we have 
no reflection. The condition aad = —bad is reached for 
Ed ~ 5.2, but the effect is even more impressive at some- 
what smaller Ed (in Fig. [1] we deliberately show the case 
Ed = 4) . Then the region of collision energies where the 
transmission probability is smaller than a certain small 
but finite value is wider. 

Finally, let us comment on the dimer-dimer collisions. 
The scattering length for two quasi- ID dimers consisting 
of fermionic atoms has been calculated in Ref. [Sg]. In 
principle, one can also find the dimer-dimer phase shift 
at finite collision energies. However, below the break- 
up threshold the two-dimcr collisions do not lead to mo- 
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mcntum diffusion even if the scattering phase shift differs 
significantly from the zeroth order one. This is because 
these are identical particles of the same mass. The mo- 
mentum diffusion should appear in three-dimer collisions 
or in other processes involving more particles (atoms 
and/or dimers). Of course, these are much more difhcult 
to analyze, but we can neglect these few-dimcr diffusion 
channels compared to the two-body atom-dimer channel 
when the density of dimers is small or when they are in 
the Tonks regime, i.e. when the probability to find two 
dimers close together is suppressed. 



VIII. SUMMARY AND CONCLUSION 

We have developed a perturbation theory for the quasi- 
ID three-atom problem in the near-integrable limit. We 
have shown that to the leading order the integrability of 
the f-t-i fermionic system is broken by the effective range 
corrections to the two-body coupling constant originat- 
ing from the virtual transverse excitations. In contrast, 
the quasi-lD problem of three bosons can be reduced 
to the purely ID problem by adding an additional local 



three-body term [11[, the two-body effective range cor- 
rections and other effects being of higher order. We have 
calculated the atom-dimer reflection, transmission, and 
break-up probabilities, as well as the three-body recom- 
bination rate constants for fermions and bosons as func- 
tions of the energy and the interatomic scattering length. 
At a finite collision energy q^ the reflection probability 
is proportional to aF' /l}^^q^ in the fermionic case and to 
a^ /l^q^ for bosons. 

We have shown that for fermions the rates of relax- 
ation and recombination to deep molecular states are 
much lower than the rates of (momentum) diffusion in 



the zeroth order Bcthe ansatz basis, meaning that the in- 
tegrability breaking processes occur well within the life- 
time of the system. In this respect the fermionic case 
differs strongly from the bosonic one. For bosons the 
rate of inelastic relaxation/recombination is comparable 
or higher than the rates of momentum diffusion. How- 
ever, the lifetime of the bosonic gas is sufficiently long 
to study "zeroth order" effects of integrability, such as 
the reflectionless atom-dimer scattering. In spite of the 
absence of reflection the atom-dimer interaction is strong 
as can be seen from the phase of the transmission ampli- 
tude. The corresponding shift in position or in time of 
the relative atom-dimer trajectory can be measured, for 
example, in a quantum Newton's cradle type of experi- 
ment pTj . 

The strong a-dependence of the atom-dimer reflection 
probability and the suppression of inelastic processes in 
a mixture of quasi-lD f-fermions and 'f-^-dimers makes 
this system an ideal candidate for a controllable investi- 
gation of differences in integrable and nonintcgrable dy- 
namics. We have estimated parameters of a realistic ex- 
periment, which can be performed in a spin-imbalanced 
mixture of two hyperfine states of ''°K. 
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